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I. INTRODUCTION 

Elucidation of acoustical mechanism of air-reed instru- 
ments is a long standing problem in the field of musical 
acoustics and is still not understood completely 1,2 . The 
sound source of those instruments is a sort of aerody- 
namic sound source, so-called edge tone, which is gener- 
ated by an oscillating jet flow collided with an edge 3,4 . 
Actually, the oscillating jet flow, which emanates from a 
flue, passes through an open mouth and collides with the 
edge, plays the role of sound source. The major difficulty 
of analysing the acoustic mechanism of air-reed instru- 
ments is in strong and complex interactions between the 
air flow dynamics acting as the sound source and a sound 
field excited in a resonator by \\^£&. Indeed, the reso- 
nance of a pipe induces the sound field exceeding 140dB 
in it so that the strong sound field must affect the motion 
of the sound source jet, if the synchronisation between 
the jet oscillation and the sound field is well sustained. In 
order to elucidate this process in detail, we need to study 
the dynamics of the jet flow in terms of fluid dynamics 
as well as that of the sound field in terms of acoustics 
at the same time taking into account the complex inter- 
action between them. However the method necessary to 
attack this problem has not been established yet. 

In the long history of studying the edge tone, some 
phenomenological formulae that describe the relation of 
the edge tone frequency with the jet velocity have been 
proposed so far. A very useful formula was introduced 
based on the experimental results by Brown in 1937 3 : 
it indicates that the frequency is proportional to the jet 
velocity. In the field of the musical acoustic, the phe- 
nomenological theory that describes the behavior of in- 
strument driven by an air jet has been developed since 
1960's 1 . The mechanism that the air jet drives the res- 
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onator was studied by several authors experimentally 
and theoretically, chief among pioneers being Cremer and 
Ising, and Cohmar^^^^^^dldld^dld^d^. It was found 
that two type of driving mechanism, volume-flow mecha- 
nism and momentum-drive mechanism, simultaneously 
work in jet-resonator interaction, though the volume- 
flow mechanism dominates in a lower range of the jet 
velocityi*£i£*iM&*i£. The theory that describes the be- 
havior of the resonator driven by the air jet has also de- 
veloped with help of the equivalent electrical circuit the- 
ory. As a result, the difference in sounding between the 
pure edge tone and the air-reed instrument is clarified in 
relation of the sound frequency to the jet velocity. 

However, the theory developed in the field of musical 
acoustics includes many conceptual approximations and 
is far from rigorous. A rigorous theory must be framed 
based upon the Navier-Stokes equations of fluid dynam- 
ics considering the role of vorticity field of a real fluid as a 
sound source in the complex geometry of the instrument. 
The aerodynamic sound source was first formulated by 
Lighthill in 1952 17 . Lighthill introduced an inhomoge- 
neous wave equation whose inhomogeneous term behaves 
as a quadrupole source term, so-called Lighthill's acous- 
tic analogy. Powell and Howe followed Lighthill's work 
and introduced a very important notion that the ma- 
jor part of Lighthill's source term comes from unsteady 
motion of vortices, namely vortex sound 1 ^^. Further, 
Howe discussed the acoustic mechanism of air-reed in- 
struments in terms of the vortex sound theory^i^, which 
has since been followed in some detail by Hirschberg et 
a j2 4 20 i 2i i 22 i 23 i 24 1 25 > Nevertheless, the detail mechanism of 
sound generation by the air jet and of jet-resonator inter- 
action is still an unsolved problem due to complex behav- 
ior of fluid in the complex geometry of the instrument, 
though great efforts have been made on this problem and 
are continued^^^^^^L 3 ^ 3 -. 

The final goal of our study is to analyze the inter- 
action of the jet flow with the sound field excited in the 
resonator and to explain in terms of the vortex sound the- 
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ory the acoustic mechanism which excites and sustains a 
sound field in the resonator of air-reed instruments 4 . For 
the first step, we numerically study sounding vibrations 
of a small air-reed instrument by using a two-dimensional 
model in this paper. Comparing with the theory devel- 
oped by Cremer & Ising, Coltman and other authors, we 
will discuss to what extent our model reproduces char- 
acteristic features of sound vibrations of air-reed instru- 
ments. 

For the numerical study on aero-acoustic problems, 
e.g., noises caused by aircraft and high-speed trains, 
we normally use a hybrid method, in which the sound 
field can be separated from the turbulence flow acting 
as sound sources and they are calculated by different 
schemes, each of which is suitable for the calculation of 
individual dynamics. However, it is necessary for the 
calculation of air-reed instruments to solve the turbu- 
lence flow and the sound field simultaneously, because 
of strong nonlinear interaction between the flow and the 
sound. Then, we use a compressible LES (Large-Eddy 
Simulation) instead of the hybrid method 34 . The reason 
we adopt LES is that LES is very stable for a long term 
calculation, though it somewhat sacrifices accuracy. 

The organization of the present paper is as follows. 
In section 2, we review previous studies relative to the 
sound production of air-reed instruments. We briefly ex- 
plain Ligthhill's theory and Brown's work on the edge 
tone. Further, we explain the phenomenological theory of 
air-reed instruments developed by Cremer & Ising, Colt- 
man and other researches. In section 3, we introduce 
our model instrument, a small two-dimensional air-reed 
instrument with a closed end, and explain the environ- 
ment of numerical calculation. In section 4, we show the 
results of the numerical analysis. First we explain re- 
sults of stable oscillation with an optimal choice of jet 
velocity. The relation of time evolution of acoustic pres- 
sure observed in the pipe with that of vorticity of the 
jet flow is investigated taking the correlation between 
them. We also show special distributions of character- 
istic dynamical- variables observed at a certain time in a 
near field of the instrument: air density, velocity, vortic- 
ity and Lighthill's sound source are considered. Next, we 
study frequency change of acoustic wave excited in the 
pipe with increase of the jet velocity. We compare our 
numerical results with the prediction given by the phe- 
nomenological theory of air-reed instruments introduced 
in section 2. Section 5 is devoted to a summary and 
discussion. 



ity equations, to an inhomogeneous wave equation whose 
inhomogeneous term plays the role of the source: 
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where the tensor Tij is called Lighthill's tensor and is 
defined by 

T^ = pViVj + (0 - p ) - c 2 (p- po))Sij + (Tij. (2) 

Here, c denotes the speed of sound, p the air pressure 
with the average po> P the air density with the average 
po, and (Tij the viscous stress tensor. The sound gener- 
ated by the localized turbulence is regarded as that wave 
propagating in a stationary acoustic medium, which is 
generated by the quadrupole source distribution given 
by the inhomogeneous term in RHS of eq.(pQ). This is 
Lighthill's acoustic analogy. 

Since the dissipation by cr^ can be ignored for a 
high Reynolds number and adiabaticity is well held as 
(p — Po) — c<2 (p — Po) = 0, then the first term of eq.(J2]), 
pViVj, becomes the major term of the source. Further, 
particle velocities of the sound are usually sufficiently 
small comparing with those of the real flow and so the 
source term is well approximated by that obtained from 
incompressible fluid with p = po and div v = 0. For two 
dimensional (2D) fluid, it is reduced into 



d 2 T tJ 
dxidxj 



-2p 



dvi dv2 dv2 dvi 
dxi dx2 dxi dx2 



(3) 



In calculation of Lighthill's source in section IIV1 we will 
use this formula. 



B. Edge tone 

As shown in FigJTJ edge tone is an aerodynamic sound 
generated by the unsteady but mostly periodical oscilla- 
tion of the jet emanated from the flue and collied with 
the edge. The edge tone is the sound source of air-reed 
instruments^^. Although the detail mechanism of the 
production of edge tone has not been completely under- 
stood yet, its characteristic features have been well cap- 
tured by semi-empirical equations introduced based on 
experimental results. One of those equations was intro- 
duced by Brown to predict the frequency of edge tone^: 



v = 0.466j(100V - 40)(1/(100Z) - 0.07), 



(4) 



II. THEORY RELATIVE TO SOUND PRODUCTION OF 
AIR-REED INSTRUMENTS 

A. Lighthill's Theory 

The sound generated by turbulence is usually called 
aerodynamic sound, which is a very small byproduct of 
the motion of unsteady flows of high Reynolds number. 
The source of aerodynamic sound was first given the ex- 
act form by Lighthill 17 . Lighthill transformed exactly the 
set of fundamental equations, Navier- Stokes and continu- 



where V denotes the speed of jet and I is the distance 
between the flue and the edge. The number j is taken as 
j = 1.0,2.3,3.8,5.4. For j = 1, it gives the fundamental 
frequency and others denote overtones. With increase 
of V, the fundamental oscillation is excited and its fre- 
quency increases in proportion to V. But it jumps to 
one of overtones, if V exceeds a threshold value, and it 
jumps successively from one to other with increase of V. 
The transitions are hysteretic so that the threshold val- 
ues of V of the downward process are usually different 
from those of the upward. 
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FIG. 1. Edge tone 



FIG. 2. Jet oscillation 



C. Phenomenological theory of air-reed instruments 

The driving mechanism of air-reed instruments was 
first studied by Cremer & Ising, and Coltman as pio- 
neer works£&£&£, which have been followed by many 
author o 14 i 15 i 16 . Those studies have made differences in 
driving mechanism clear between the edge tone and the 
air-reed instrument with a resonator. It is turned out 
that there are two ways in which a pulsating jet drives 
a pipe, momentum drive and volume flow drivei. In the 
momentum drive, the jet is brought to rest by mixing 
and dissipation, and thereby generates an acoustic pres- 
sure which acts on an acoustic flow in the resonator. In 
the volume flow drive, the jet essentially contributes a 
volume flow which acts on the acoustic pressure at the 
open mouth, which is not zero because of the end correc- 
tion. In the following, we mention the outline of driving 
mechanism and regenerative excitation mechanism of air- 
reed instruments according to the text book by Fletcher 
and Rossing 1 . 



1. Driving mechanism 

The jet which emanates from a flue is disturbed by 
a uniform transverse acoustic-flow v z ex.-p(iujt) and forms 
an oscillating wave which propagates with a velocity u 
and grows exponentially with a growth parameter fi. The 
wave of jet propagating in x-direction is well approxi- 
mated by the semi-empirical equation (also see Fig]2]): 

J(x) = — i — {ex.p(iujt) — cosh(fix) ex.p(iuj(t — x/u))}. 

UJ 

(5) 

In the range kd ^> 1, where d stands for the height of flue 
aperture, we can use approximations u ~ V/2 and \i ~ k. 
Since they are available in a quite wide range, we always 
assume those approximations in the following. In this 
formula, the existence of edge, which plays the important 
role in formation of the jet oscillation, is ignored. Thus 
it should be considered that eq.flSJ) only gives the lowest 
order approximation of the jet motion. 

Taking into account the jet oscillation given by eq.flSJ), 
it is possible to make phenomenological formulation of 
that driving mechanism through which the jet acts on 
the pipe. As shown by Figj3j after the jet of a steady 
velocity V reaches the edge of the pipe, a fraction of its 
cross section Sj enters the pipe. It is assumed that the 



fraction of jet blends with an acoustic flow in the volume 
between the cross sections M and P with a distance Ax 
so that only acoustic motions survive in the other part of 
pipe beyond the cross section P. Under this assumption, 
the acoustic volume flow U p at the cross section P is 
expressed by 



u„ 



pV 2 Sj 



ipujALVSj 
Sp[Z p + Z m ) Dp(Z p + Z m ) 



(6) 



where Z m = ipooAL/ 'S p denotes the impedance at a fre- 
quency uj looking out from the plane M, AL the end 
correction at the open mouth, Z p the pipe impedance 
evaluated at the plane M and Sj the area of the cross 
section of the pipe. 

Two term in RHS of eq. (|6j) indicate two different drive 
mechanisms. The first constitutes the momentum drive 
and the second is the volume flow drive. In the situation 
ujAL > V, the second term is larger than the first so that 
the volume-flow mechanism dominates the momentum- 
drive mechanism. For our numerical calculations imple- 
mented in sectiorflVl substitution of representative values 
of the fundamental frequency of the pipe and the end cor- 
rection of the open mouth, i.e., uj ~ 2ir x 913.1rad/s and 
AL = 0.005m, gives ujAL = 9.131tt - 28.7m/s. Thus 
the volume-flow mechanism dominates if V < 28.7m/s. 
For the third harmonic at 2739.3Hz, the volume-flow 
mechanism governs the driving process in the range 
V < 86.1m/s. As shown later, the fundamental domi- 
nates for V < 22m/s, while the third harmonic overcomes 
the fundamental for V > 24m/s, then it is expected that 
the dominant mechanism exciting the sound wave in our 
numerical range (2 < V < 40m/s) is the volume flow 
drive. 
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FIG. 3. Pipe driven by jet 
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FIG. 4. Equivalent circuit network 



2. Regenerative excitation mechanism 

Based on the analysis in the previous subsection, we 
can describe the sound excitation mechanism of air-reed 
instruments. To do this, we separate the system into a 
generator and a resonator by the plane, such as P or M. 
Make use of an equivalent circuit network (see FigH]), 
we can determine the oscillating condition giving sound 
frequencies of the modes as functions of the jet velocity. 

The stability condition of the network is given by 



Yj 



* m i J- n 



o, 



(7) 



where Yj denotes the jet admittance obtained from the 
jet motion in eq.©, which is given by 



Yi 



VW 
ptu 2 AL 



cosh jil exp 



(?+*) 



(8) 



satisfied, we obtain kL ~ nir for the open end pipe and 
kL ~ (n + l/2)7r for the closed end pipe, which give the 
resonance conditions for the open and closed end pipes, 
respectively. 

Fig|5] illustrates the frequency of excited sound wave 
as a function of the jet velocity V given by eq. (fTT|) for 
the closed end pipe together with the ideal condition 
ujI/u = 7r — (j) and the edge tone frequency in eq.fllj) with 
j = 1. The frequency of sound wave first increases in 
proportion to the jet velocity similar to the edge tone, 
but it is synchronized with the fundamental resonance, 
if the edge tone frequency comes near the fundamental 
frequency. Locking to the fundamental is continued until 
the edge tone frequency is close to the next resonance 
frequency. After that it jumps up to the next resonance 
and is synchronized with it, so the same process may 
be repeated every when the edge tone frequency reaches 
next one. 



-*^ 




Jet Velocity 



with (j) = arctanl ^7; ) ? the mouth admittance Y m is de- 
fined by Y m = 1/Z m = -z T L , and the pipe admittance 

Y p has the forms: Y p = — i-£ cot kL for an open end pipe 

s 
with a length L and Y p = z— tan&L for a closed end 

pipe. 

It is convenient to separate the condition ([7]) into the 

real and imaginary parts 



ReYj < 0, 



TmfYj + Y m + Y p 



0. 



(9) 
(10) 



Eq.(|9j) indicates that the jet must have a negative re- 
sistance as a power supplier that excites the system 
overcoming non-zero resistances of the mouth and pipe 
impedance, which are, however, ignored in this analy- 
sis. For the closed end pipe, the imaginary part (|TD|) is 
rewritten as 



1 



V 



jAL 



j 2 ALh 



cosh/i/ sin 



sujI ( V 
h arctan - 



\ujAL 



-tanfcL = 0. (11) 

c 



The ideal excitation is given by ReK, < and ImYj = 
so that the phase of Yj in eq.© takes — 7r, i.e., ujI/u = 
7r — (j) and Im(Y m + Y p ) — 0. Since ImF m < 0, Iml^, must 
takes a positive value. If the condition uoAL/c <C 1 is 



FIG. 5. Oscillation frequency vs. jet velocity 



III. MODEL AND NUMERICAL SCHEME 

A. Numerical model 

For numerical analysis of air-reed instruments, we si- 
multaneously calculate the dynamics of the jet flow and 
the sound field excited in a resonator. The sound speed c 
of about 340m/s is much higher than the jet velocity V, 
which is at most several tens in MKS units. For reproduc- 
tion of sound, an extremely smaller time step is required 
compared with ordinary numerical calculations of fluid 
dynamics. On the other hand, spatial scales used for cal- 
culations of fluid dynamics with those vortices, some of 
which may be smaller than 1mm, are much less than wave 
lengths of sound, e.g., 34mm even at 10kHz. Therefore, 
in the numerical calculation of air-reed instruments we 
must satisfy the both requirements, a sufficiently small 
time step to describe sound propagation and a spatial 
mesh finer enough to reproduce vortices in fluid. Fur- 
ther, particle velocities of sound (or energies of sound) 
are usually much less than those of the flow. It turns out 
that sound energies in living environment are 10 -4 times 
as small as or smaller than those of fluid. It is not easy 



Numerical study on an air-reed instrument with LES 4 



(a) 



EL 



d = l 



-&L 



^L 



6=25 



1 = 5 



h = 10 * 



L=90 



(b) 



i I 



FIG. 6. (color online) Model and mesh, (a) Dimensions of 2D model (Unit of length :mm; Unit of angle: degree), (b) Numerical 
mesh. Dotted lines denote transparent walls. 



TABLE I. Parameters of mesh 



points 


cells 


faces 


158,762 


78,492 


314,856 



to numerically calculate with a high degree of accuracy 
sound propagation dissipating for a long distance. 

To realize the calculation, we take a two dimensional 
(2D) model of a small air-reed instrument and concen- 
trate our attention on dynamics in a near field of it. Tak- 
ing a 2D model rather than a 3D model makes the num- 
ber of grids necessary for the calculation reduce, which 
gives rise to diminishing elapse time, although some of 
important 3D effects might be lost. 

Fig ED (a) shows geometry of the 2D model we adopt, 
where d denotes nozzle height, I width of mouth aper- 
ture, L pipe length, h pipe height and edge angle, and 
those values are taken as shown in the figure. The model 
instrument is 9mm in length thereby the fundamental 
frequency being estimated as Jo °° 913.1Hz including the 
open end correction. The edge angle is fixed at 25°, at 
which we got the most stable oscillation in preliminary 
calculations. This angle value is in a suitable range for 
real instruments. 



B. Numerical method 

For numerical calculation, we use a compressible LES 
(Large eddy simulation), which is very popular in numeri- 
cal simulations of aero-acoustics. Actually the scheme we 
adopt is the compressible LES solver, Coodles of Open- 
FOAM, provided as a free software by OpenCFD Ltd 35 . 
LES is very stable for a long time simulation, while it 
involves some ambiguities in boundary layers due to sta- 
tistical assumption for dynamics of eddies smaller than a 
given grid size. It however makes sense in the statistical 
point of view. 

As shown in FigEfb), the spatial area of the numer- 
ical mesh is taken large enough to include an acoustic 
field outside the instrument and parameters of the mesh 
are taken as shown in Table 1. The averages of pres- 



sure and temperature are taken as po = lOOkPa and 
To = 300K, respectively. The time step of numerical in- 
tegration is At = 10 _7 sec and time evolution up to O.lsec 
is obtained. The velocity of the jet emanating from the 
flue V is changed as a control parameter in the range 
(2 < V < 40m/s). Observations are done as follows. 
The acoustic pressure p and air density p are observed at 
the point A, distance of 10mm from the right end of the 
pipe and on the center line of it. The vorticity uo in the 
jet flow is calculated at point B, distance of 1.6mm right 
from the exit of the flue and on the center line of it. 



IV. NUMERICAL RESULTS 
A. Stable oscillation 

In this subsection, we show numerical results of that 
oscillation at V = 12m/s which is the most stable among 
oscillations observed in the whole range of V. Figd (a) 
and (b) show the acoustic pressure oscillation observed at 
the point A and its power spectrum, respectively. For the 
calculation of the spectrum, initial transient oscillations 
(0 < t < 0.01s) are omitted. The oscillation of the acous- 
tic pressure in FigJT] (a) is very stable with an almost 
constant pitch in the whole range of time (0 < t < 0.1s) 
except for a short initial transition. Its amplitude is gen- 
tly undulated but reaches up several hundreds Pa, which 
are much larger than normal acoustic pressures in the 
open air field. As shown in FigJTJb), the main peak of 
the spectrum is very sharp and appears at v = 805.7Hz, 
which is less than the theoretical estimation of the funda- 
mental resonance frequency 913.1Hz by 107.4Hz. How- 
ever, it rapidly approaches the resonance frequency with 
increase of V as shown later. Thus it is considered that 
the oscillation at V = 12m/s is almost locked on the 
fundamental resonance. 

FigH] (a) and (b) show the vorticity fluctuation of jet 
observed at the point B and its power spectrum, respec- 
tively. The vorticity regularly oscillates like the acoustic 
pressure at the point A with the same fundamental fre- 
quency as shown in the power spectrum(see Fig {8] (b)). 
Fig |9] shows the mutual correlation function between the 
acoustic pressure (0.01 < t < 0.1s) and the vorticity 
(0.01 < t < 0.1s). The correlation changes periodically 
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FIG. 7. Acoustic pressure at point A. (a) Acoustic pressure 
oscillation, (b) Power spectrum. 



FIG. 8. Vorticity at point B. 
Power spectrum. 



(a) Vorticity oscillation, (b) 



with the same fundamental frequency and does not decay 
for a long time. It means that there is a strong correlation 
between the jet motion and sound field in the resonator 
and so the jet motion is controlled rather by the pipe 
resonance than by its self-oscillation mechanism. 

FiglTUl shows spacial distributions of characteristic dy- 
namical variables at a certain time: air density, flow ve- 
locity, vorticity and Lighthill's sound source. We used 
eq.(|3J) for the calculation of Lighthill's sound source. 
Though it is not shown in a picture, the acoustic pressure 
distribution p(x, t) is almost as same as the air density 
distribution p(x,t). In the stationary oscillation regime, 
fluctuations of p(x,t), more precisely \p(x,i) — po\, be- 
come much larger in the inside of pipe than in the outside 
as shown in FigJlO] (a). 

In the velocity distribution in Fig lToT b), we see that 
eddies mainly arise from collision of the oscillating jet 
with the edge and are pouring into the inside or going to 
the outside of pipe. The eddies in the outside run along 
the wall of pipe being dissipated gradually. On the other 
hand, those in the inside roll up and make a large rotor or 
a few rotors near the mouth opening not spreading into 
more right hand side. No eddies apparently appear in 
the right 3/4 part of the pipe, in which the strong sound 
field dominates the flow. Such a rotor (s) always appears 
in stable oscillations and so the existence of the rotor (s) 
may affect the jet motion, though we do not touch details 
of the mechanism in this paper. 

As shown in FiglToTc), vorticity takes large (absolute) 
values along the upper and lower boundaries of jet and in 
areas of strong eddies rolled up. Lighthill's sound source 
in FiglToTd) almost overlaps with the vorticity distribu- 
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FIG. 9. Correlation between sound pressure and jet vorticity. 



tion. This fact is supported by the main claim of the 
Powell-Howe vortex sound theory ^ 18 i 19 that the major 
part of Lighthill's sound sources is of contribution of vor- 
tices. The sound sources along the jet periodically os- 
cillate reflecting the motion of jet. On the other hand, 
the other sound sources mainly located on large eddies 
behave rather irregularly in time evolution due to irreg- 
ular motions of the eddies. This fact combined with 
the fact that the vorticity of jet is strongly correlated 
with the sound field excited in the pipe suggests that 
the main mechanism driving the sound field in the pipe 
at V = 12m/s is 'volume flow drive' rather than 'mo- 
mentum drive'. This suggestion is also supported by the 
theoretical prediction discussed in subsection lILC.ll 
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FIG. 10. (color online) Spatial distributions of representative dynamical variables, (a) Air density, (b) Flow velocity, (c) 
Vorticity. (d) LighthiU's sound source. 



B. Frequency change with jet velocity 

In this subsection, we discuss change of characteris- 
tic frequencies of sound waves excited in the pipe, i.e., 
spectrum peaks of fundamental and first overtone, with 
increase of the jet velocity V. In FigJTTJ we show charac- 
teristic frequencies of the acoustic pressure at the point 
A as functions of V in the range (2 < V < 40m/s). For 
comparison, resonance frequencies of the pipe estimated 
theoretically and the edge tone frequency given by eq. (|4]) 
with j = 1 are also depicted in FigJTTJ The fundamen- 
tal frequency of the acoustic pressure is observed in the 
whole range of V, but the first overtone is identified only 
at high velocities, V > 18m/s. 
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FIG. 11. (color online) Changes of oscillation frequencies with 
jet velocity. The lines labeled 'pipe res.l' and 'pipe res. 3' 
respectively indicate the frequencies of fundamental and 3rd 
harmonics resonances estimated theoretically. 



In the low velocity regime, (V < 8m/s), the fundamen- 
tal frequency of acoustic pressure increases as that of the 
edge tone, namely it is proportional to V. Hence the jet 
motion is little affected by the pipe resonance and almost 



keeps its natural oscillation. For example, the acoustic 
pressure and its power spectrum at V = 4m/s are shown 
in FiglT^l (a) and (b), respectively. The acoustic pressure 
oscillates somewhat regularly at vq = 341.8Hz with very 
small amplitudes compared with those at V = 12m/s. 
However, the power spectrum shown in Figjl2^b) is some- 
thing interesting: the first overtone peak at v\ = 659.2Hz 
accompanied by smaller side peaks in a upper range is ob- 
served. This broad band of peaks seems to be affected by 
the pipe resonance. Indeed, the fundamental frequency 
vq is slightly less than half of fundamental resonance 
frequency, so a broad band of peaks accompanying the 
overtone can appear near the fundamental resonance fre- 
quency. As shown in FiglT^l (c). the correlation between 
the acoustic pressure and vorticity oscillates in a regu- 
lar manner but gradually decays to an oscillation much 
smaller than that at V = 12m/s. Then it is considered 
that the interaction between the generator and resonator 
is small in this regime. Namely an unidirectional inter- 
action from the generator to the resonator dominates the 
operation with less feedback from the resonator, thereby 
being out of resonance. 

In the middle range, (10 < V < 22m/s), oscillations of 
sound wave locking on the fundamental pipe resonance 
are observed. The frequency of oscillation is fairly lower 
than that of fundamental resonance at V = lOm/s, but it 
quickly approaches the frequency of the resonance with 
increase of V. The oscillations are very stable in the 
range (10 < V < 16m/s) and the most stable one is 
observed at V = 12m/s as shown in the previous subsec- 
tion. So the oscillation is most stabilized just above the 
value of V at which the jet motion starts synchronizing 
with the pipe resonance. In the range (18 < V < 22m/s), 
oscillations, however, become slightly unstable. Namely 
some amplitude modulations occur in long term evolu- 
tions, though we do not show results. In the spectra of 
those oscillations, there observes a peak corresponding to 
the edge tone frequency at the present velocity V. Fur- 
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FIG. 12. Oscillation at V = 4m/s. (a) Acoustic pressure, (b) FIG. 13 
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Oscillation at V = 36m/s. (a) Acoustic pressure, (b) 
(c) Correlation with vorticity at point B. 



ther smaller peaks whose frequencies are nearly equal to 
the third harmonic of the pipe are often observed. It is 
considered that competition among the resonances and 
edge tone, i.e., inherent oscillation of the jet, occurs, but 
the fundamental still dominates the others, though its 
oscillation is somewhat disturbed. 

In the high velocity range, (V > 24m/s), first over- 
tone peaks likely to be the third harmonic are clearly 
observed in the spectra of acoustic pressure. The fre- 
quency of the first overtone peak increases mostly pro- 
portional to the jet velocity like the edge tone in the 
range 18 < V < 24m/s but almost converges to a con- 
stant value in the range V > 32m/s, e.g., v\ = 2490.2Hz 
at V = 36m/s, although it is considerably lower than the 
theoretical estimation of the third harmonic, 2739.3Hz. 
Then the oscillation would start to be synchronized with 
the third harmonic around V = 24 ~ 26m/s. Indeed, 
the peak height of the first overtone is increasing with 
increase of V and becomes larger than that of the fun- 
damental for V > 24m/s (for example see Fig fT3T b) at 
V = 36m/s). Further the wave form of acoustic pressure 
changes to that of the third harmonic for V > 24m/s. 



Fig IT3t a) shows the time evolution of the acoustic pres- 
sure at V — 36.0m/s. It oscillates almost periodically 
with very strong amplitudes, sometimes exceeding lkPa. 
As shown in FigHHTc), the correlation between the sound 
pressure at the point A and the vorticity at the point B 
behaves in a regular manner and does not decay in a 
long term. In conclusion, it can be said that the transi- 
tion from the fundamental to the third harmonic occurs 
around V = 24m/s. 

We compare our numerical result with the theoretical 
prediction introduced in section Hi. C. 2 1 FigIT4l shows the 
velocity- frequency curves (broken lines) given by eq. ([TT]) 
at parameter values adjusted to our numerical calcula- 
tion: I = 5mm, AL = 9mm, h = 10mm, L = 90mm, 
c = 340m/s, fil = 7r. However, we need to rescale the ar- 
gument of sinusoidal function in eq. ([TT]) as u/l = 0.8V/ 1 
to make it be adjusted to the edge tone equation (|4]) and 
also to our numerical result, while it is usually taken as 
u/l < 0.5V '/I. The curves of eq. ([TTJ ) are drawn only in the 
area limited by the two dotted lines given by Relj = 0, in 
which the jet has a negative value in resistance, the nec- 
essary condition to excite the instrument. The chain line 
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FIG. 14. Comparison of the numerical result with the theoret- 
ical prediction in change of frequencies with the jet velocity. 



stands for the optimum oscillation condition Imlj = 0. 
For comparison, the velocity-frequency curves obtained 
numerically, which are the sames as those in FigJTTJ are 
also drawn by solid lines. 

Our numerical result shows good agreement with the 
theoretical prediction in three characteristic ranges: os- 
cillations of the edge tone in the low velocity range, lock- 
ing to the fundamental resonance in the middle range, 
and transition and locking to the third harmonic in the 
high range. It is expected that the optimum oscillations 
of the first and third harmonics occur at intersections 
of the corresponding branches of eq. (fTT|) with the line 
Imlj = 0, respectively. The intersection appears just af- 
ter locking to the first or third harmonics starts. The 
fact that the optimum oscillation of the fundamental is 
observed numerically at V = 12m/s well supports the 
validity of theoretical prediction. A well sustained third 
harmonic, e.g., at V = 36m/s, is also observed near the 
intersection of the third harmonic. 

Here, we explain why we take the parameter u/l as 
u/l = 0.8V// instead of u/l = 0.5V// in order to cal- 
culate the velocity- frequency curves and optimum line. 
According to experiments and the semi-empirical theory, 
the velocity of jet wave u should take a value in the range 
u < 0.5V. However, even at u = 0.5V, the velocity- 
frequency curves given by eq. ([TT]) lean to the right more. 
As a result, the curves are markedly shifted from the edge 
tone line in a low velocity range and the onset point of 
locking to the fundamental is estimated as V ^ 20m/s. 
Then those curves apparently go out of our numerical 
results and also out of Brown's edge tone equation. Tak- 
ing into account the fact that the existence of the edge 
is ignored in framing the jet model in that theory (see 
IIIC.l|h it is considered that such extent of modification 
is still within an acceptable limit. It is rather surprised 
that the semi-empirical theory shows quantitatively good 
agreement with full numerical calculations with such a 
slight modification. 

Finally note that as discussed in subsection IILC.1| it 
is predicted that the volume-flow mechanics dominates 
in excitation of sound waves not only in the range of the 
fundamental but also in the range of the third harmonic. 
We think that this prediction is probably true and is 
partially supported by the observation that the correla- 



tion between the acoustic pressure and the vorticity of 
the jet behaves regularly without remarkable decays in 
the ranges of V, in which the fundamental and the third 
harmonic are well sustained. Further it is also confirmed 
in preliminary calculations that the correlation becomes 
more unstable, if the vorticity is observed at a point in 
the relaxation area between the planes M and P in Fig]3j 
i.e., momentum driving area, because motions of eddies 
in this region are more irregular than the jet oscillation. 



V. SUMMARY AND DISCUSSION 

In this paper, we have reported on the numerical anal- 
ysis of a 2D air-reed instrument with the compressible 
LES. As a result, vibrations of the air-reed instrument 
are well reproduced by using the compressible LES. Espe- 
cially, some characteristic features of the air reed instru- 
ment obtained numerically show good agreement with 
those given by the theoretical prediction as well as experi- 
mental results^&£ i 14 i 15 i 22 . The characteristic features re- 
produced are as follows. 

When the velocity is small enough, the inherent jet os- 
cillation radiating the edge tone is predominant so that 
the oscillation frequency of sound excited in the pipe is al- 
most proportional to the jet velocity. However, the jet os- 
cillation is synchronized with the fundamental resonance 
of the pipe when the frequency of edge tone approaches 
that of the fundamental resonance with increase of the jet 
velocity. The synchronization with the fundamental res- 
onance continues until the frequency of edge tone reaches 
that of the first overtone, i.e., 3rd harmonic in our case. 
Then the synchronization is turned off and the transition 
to the 3rd harmonic arises, namely the frequency locking 
to the 3rd harmonic takes place. 

Based on the results of this paper, we are able to pro- 
ceed to study the acoustic mechanism of the air-reed in- 
strument comprehensively. With help of the Lighthill 
theory 17 and/or Powell-Howe vortex sound theor y 19 ! 18 ^, 
the place of sound sources in turbulence will be detected 
clearly and their behavior will be characterised in terms 
of aero-dynamics and nonlinear dynamics. The follow- 
ing problems and questions should be investigated. The 
questions which mechanism, volume-flow mechanism or 
momentum drive mechanisn>^^^ i 14 i 15 i 16 , dominates at 
a given playing condition and what change occurs with 
change of the playing condition, e.g., with the jet veloc- 
ity, should be answered clearly. It should be also clarified 
what kind of influence the resonance of the pipe exerts 
on the motion of the jet, because it must be the key to 
understand the synchronization mechanism between the 
jet flow and pipe resonance, the key mechanism of wind 
instruments different from the edge tone. 

To analyze more realistic behavior of instruments, the 
numerical analysis of 3D models is needed. In a prelim- 
inary calculation, we have succeeded to numerically re- 
produce acoustic oscillations of the ocarina with a 3D 
model 36 . It is considered that the ocarina uses the 
Helmholtz resonance caused by an elastic property of 
air instead of the pipe resonance. Then, it is quite in- 
teresting to consider the problem how different acoustic 
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mechanisms work for different types of resonator with 
comparing numerical data of the ocarina with that of a 
3D model of the instrument with a resonance pipe. Be- 
fore that it is of course necessary to clarify differences 
between the 2D and 3D models. 
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